Free stream capturing in fluid conservation law for moving coordinates in three dimensions

The free-stream capturing technique for both the finite-volume (FV) and finite-difference (FD) framework is summarized. For an arbitrary motion of the grid, the FV analysis shows that volumes swept by all six surfaces of the cell have to be computed correctly. This means that the free-stream capturing time-metric terms should be calculated not only from a surface vector of a cell at a single time level, but also from a volume swept by the cell surface in space and time. The FV free-stream capturing formulation is applicable to the FD formulation by proper translation from an FV cell to an FD mesh

Use of high-resolution upwind scheme for vortical flow simulations

For vortical flow simulations at high Reynolds numbers, it is important to keep the artificial dissipation as small as possible since it induces unphysical decay of the vortex strength. One way to accomplish this is to decrease the grid spacing. Another way is to use computational schemes having little dissipation. Here, one of the high-resolution upwind schemes called MUSCL with Roe's average is applied to vortical flow fields. Two examples are considered. One is the leading-edge separation-vortex flow over a strake-delta wing. The other is a high-angle of attack supersonic flow over a spaceplane-like geometry. Comparison with the central difference solutions indicates that the present upwind scheme is less dissipative and thus has better resolution for the vortical flows

Turbulent jet interaction with a long rise-time pressure signature

A sonic boom signature with a long rise time has the ability to reduce the sonic boom, but it does not necessarily minimize the sonic boom at the ground level because of the real atmospheric turbulence. In this study, an effect of the turbulence on a long rise-time pressure signature was experimentally investigated in a ballistic range facility. To compare the effects of the turbulence on the long and short rise-time pressure signatures, a cone-cylinder projectile that simultaneously produces these pressure signatures was designed. The pressure waves interacted with a turbulent field generated by a circular nozzle. The turbulence effects were evaluated using flow diagnostic techniques: high-speed schlieren photography, a point-diffraction interferometer, and a pressure measurement. In spite of the fact that the long and short rise-time pressure signatures simultaneously travel through the turbulent field, the turbulence effects do not give the same contribution to these overpressures. Regarding the long rise-time pressure signature, the overpressure fluctuation due to the turbulence interaction is almost uniform, and a standard deviation 1.5 times greater than that of the no-turbulence case is observed. By contrast, a short rise-time pressure signature which passed through the same turbulent field is strongly affected by the turbulence. A standard deviation increases by a factor of 14 because of the turbulence interaction. Additionally, there is a non-correlation between the overpressure fluctuations of the long and short rise-time pressure signatures. These results deduce that the length of the rise time is important to the turbulence effects such as the shock focusing/diffracting

A streamwise upwind algorithm applied to vortical flow over a delta wing

Improvements were made to a streamwise upwind algorithm so that it can be used for calculating flows with vortices. A calculation is shown of flow over a delta wing at an angle of attack. The laminar, thin layer, Navier-Stokes equations are used for the calculation. The results are compared with another upwind method, a central differencing method, and experimental data. The present method shows improvements in accuracy and convergence properties

Algorithm and code development for unsteady three-dimensional Navier-Stokes equations

A streamwise upwind algorithm for solving the unsteady 3-D Navier-Stokes equations was extended to handle the moving grid system. It is noted that the finite volume concept is essential to extend the algorithm. The resulting algorithm is conservative for any motion of the coordinate system. Two extensions to an implicit method were considered and the implicit extension that makes the algorithm computationally efficient is implemented into Ames's aeroelasticity code, ENSAERO. The new flow solver has been validated through the solution of test problems. Test cases include three-dimensional problems with fixed and moving grids. The first test case shown is an unsteady viscous flow over an F-5 wing, while the second test considers the motion of the leading edge vortex as well as the motion of the shock wave for a clipped delta wing. The resulting algorithm has been implemented into ENSAERO. The upwind version leads to higher accuracy in both steady and unsteady computations than the previously used central-difference method does, while the increase in the computational time is small

Progress in computational unsteady aerodynamics

After vigorous development for over twenty years, Computational Fluid Dynamics (CFD) in the field of aerospace engineering has arrived at a turning point toward maturity. This paper discusses issues related to algorithm development for the Euler/Navier Stokes equations, code validation and recent applications of CFD for unsteady aerodynamics. Algorithm development is a fundamental element for a good CFD program. Code validation tries to bridge the reliability gap between CFD and experiment. Many of the recent applications also take a multidisciplinary approach, which is a future trend for CFD applications. As computers become more affordable, CFD is expected to be a better scientific and engineering tool

Algorithm and code development for unsteady three-dimensional Navier-Stokes equations

Aeroelastic tests require extensive cost and risk. An aeroelastic wind-tunnel experiment is an order of magnitude more expensive than a parallel experiment involving only aerodynamics. By complementing the wind-tunnel experiments with numerical simulations, the overall cost of the development of aircraft can be considerably reduced. In order to accurately compute aeroelastic phenomenon it is necessary to solve the unsteady Euler/Navier-Stokes equations simultaneously with the structural equations of motion. These equations accurately describe the flow phenomena for aeroelastic applications. At ARC a code, ENSAERO, is being developed for computing the unsteady aerodynamics and aeroelasticity of aircraft, and it solves the Euler/Navier-Stokes equations. The purpose of this cooperative agreement was to enhance ENSAERO in both algorithm and geometric capabilities. During the last five years, the algorithms of the code have been enhanced extensively by using high-resolution upwind algorithms and efficient implicit solvers. The zonal capability of the code has been extended from a one-to-one grid interface to a mismatching unsteady zonal interface. The geometric capability of the code has been extended from a single oscillating wing case to a full-span wing-body configuration with oscillating control surfaces. Each time a new capability was added, a proper validation case was simulated, and the capability of the code was demonstrated

Navier-Stokes simulation of wind-tunnel flow using LU-ADI factorization algorithm

The three dimensional Navier-Stokes solution code using the LU-ADI factorization algorithm was employed to simulate the workshop test cases of transonic flow past a wing model in a wind tunnel and in free air. The effect of the tunnel walls is well demonstrated by the present simulations. An Amdahl 1200 supercomputer having 128 Mbytes main memory was used for these computations